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ABSTRACT. A rigorous modal analysis of lamellar gratings, i.e., 
gratings having rectangular grooves, in conical mountings is presented. It is 
an extension of the analysis of Botten et aL which considered non-conical 
mountings. A key step in the extension is a decomposition of the 
electromagnetic field in the grating region into two orthogonal components. 

A computer program implementing this extended modal analysis is capable of 
dealing with plane wave diffraction by dielectric and metallic gratings with 
deep grooves, at arbitrary angles of incidence, and having arbitrary incident 
polarizations. Some numerical examples are included. 

1. INTRODUCTION 

The modal approach has been applied by many authors to lamellar, non-perfectly 
conducting gratings in the past [1-10]. Most noticeably, Botten et al. [5-7] presented a series 
of three papers, in 1981, on the modal analysis of dielectric, finitely conducting, and highly 
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conducting lamellar gratings. Their work was later formulated in a more systematic way, and 
its certain numerical aspects were improved, by Suratteau et al. [8] and by Tayeb and Petit 
[9]. However, these analyses are limited to non-conical mountings. In many applications, 
lamellar gratings are used in conical mountings [11]. Recently, a modal analysis of lamellar 
gratings in conical mountings was presented by Peng [12], 

The present work differs substantially from that of Peng [12] in mathematical 
formulation and numerical implementation. In Peng’s work, the validity of the 
decomposition of the electromagnetic field in the corrugated region into two orthogonal 
components is assumed without proof. The eigen-functions (the modal fields in the 
corrugated region) for a grating in a conical mountings are constructed by geometrical means 
from the TE and TM eigen-functions for the grating in a equivalent non-conical (also called 
classical) mounting. The completeness and orthogonality of the eigen-functions assembled 
in this manner are not addressed. In this paper, the eigen-functions along with their 
completeness and orthogonality are derived rigorously and systematically from the boundary- 
value problems. 

The present work can be considered as an extension of the works of Botten et al. and 
Suratteau et al. to conical mountings. A key step in its development is the proof of the field 
decomposition mentioned above. Once this is done, the task of finding eigenvalues and 
eigen-functions for a conical mounting reduces to that of a classical mounting, and the 
previous results of the above authors, including their powerful and sophisticated numerical 
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methods for finding the eigenvalues, can be used. 


The mathematical formulation of the modal analysis is presented in section 2, and the 
numerical aspects of the analysis are addressed in section 3. Section 4 provides some 
numerical results. The validity of the field decomposition is proved in Appendix A For the 
sake of the normal flow of the paper, some of the results available in Ref. 5-9 are re-stated 
in sections 2 and 3, but often are formulated differently. In other instances, the reader is 
referred to the original references. 

2. MATHEMATICAL FORMULATION 

2.1. Notation 

A lamellar grating in a conical diffraction configuration is depicted in Fig. 1. The 
coordinate system is chosen such that the x-axis is perpendicular to, and the z-axis is parallel 
to, the grating grooves, and the y-axis is the normal of the overall structure. A 
monochromatic plane wave of vacuum wavelength A. 0 is incident on the lamellar grating at 
a polar angle 0 and an azimuthal angle $. The range of 0 is 0 £ 0 < r/2 and that of <p 
is -ir < 0 £ ir, with the clockwise direction being the positive direction for <p. The same 
conventions will be used for the diffracted waves that are not shown. The incident 
polarization is in general elliptical. 

The geometry of the lamellar grating is shown in Fig. 2. The grating period is d, and 
the widths of medium 1 and medium 2 are d x and <f 2 * shall call the regions of space 


where y > h/2, y < - h/ 2, and - h/2 < y < h/2, regions 1, 2, and 0, respectively, where h is 

the grating groove depth. We shall use the superscript (/), where j = 1, 2, to denote 

quantities associated with regions j and the subscript j to denote quantities associated with 

the two media in region 0. Thus, the permittivity and permeability of the medium in region 

1 are c (1) and and those in region 2 are e (2) and m ( 2) - The permittivity and permeability 

of region 0 are periodic functions of x, 

€(*)=€,, li(x) - P,, 0s|x|srfj/2, 1 

e(x) - Cj, n(x) - |i 2 , dJ 2 < |x| id/ 2 . 

Although for most optical applications, the permeability is a constant and equals that of the 
vacuum, to reveal the symmetry of the electric and magnetic fields, and arc formally 
assumed to be different. The Gaussian system of units is used in this paper. 

Let k 0 be the magnitude of the vacuum wave-vector. The magnitudes of the wave- 
vector in regions j — 1, 2 and region 0 are denoted by kV and k(x), 

*» l - (2) 

t l (x) (2) 

The wave-vector of the incident plane wave is 

£ - ^(isinScoe^ - )>co®8 + t sin6 sin4> ) ■ 

We denote the z-component of the incident wave-vector by k p 

k t = * (1) sm0sm«^. ( 5 ) 

and we define the reduced magnitudes of the wave-vectors by and k(x), 
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£ U)2 - k™ - k 2 , 
P(x) = k 2 (x) - k 2 . 


( 6 ) 

(?) 


Since the physical problem is time and z-invariant, the electromagnetic field may be 
assumed to be of the form 

E(x,y,z,t) - Et,x,y)e“*~‘“‘, 

)- fay) (8b) 

Throughout this paper, i = 7-1. Substituting (8a) and (8b) into Maxwell’s equations, we 
may express the transverse components of the electromagnetic field in terms of the 
longitudinal (z) components, 


£f.x,y) - — i- [ 0 - iti) ,£x V,/f,(*oO ], 

* (X) 

fry) * *k 0 txV l E/;x,y) ], 

k (x) 


(9a) 


(9b) 


where V, = kd x + fdy. So in solving the conical diffraction problem, it is only necessary to 
work with the z-components of the electric and magnetic fields. 


As usual, the fields above and below the corrugated region may be written in 
Rayleigh expansions [13]: 

Efij) - ,<■ < IOa > 


. £ J C‘ 




(10b) 
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for y > h/2, and 


E t (x,y) - £ 




for y < /i/2, where 

a, = a 0 + 2nn/d, a 0 =* (1) sin0coe4>, 

pC/)2 ,£</)*_ a J f ^[p^] + /«[pi°] > 0. 


(11a) 

(lib) 


( 12 ) 

(13) 


In (10) and (11), I z {e) , I z (h \ Rj- e \ Rj h \ Tf \ and are, respectively, the complex 
amplitudes of the z -components of the incident and diffracted electric and magnetic fields 
in regions 1 and 2. 


2.2. Field Decomposition 

Suppose h tends to infinity in Fig. 2, yielding a medium that is periodic in x and 
infinite in y and z. If the general expression for the electromagnetic field in this infinite 
medium is known, by imposing the interface conditions aty = ± h/2 between this field and 
those given in (10) and (11), we can determine all the unknown field amplitudes. 

For a z- and y-invariant medium, the following system of equations for the fields E z 
and H z can be derived from the original Maxwell’s equations 
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w 

' t 1 d 

€ d] 


' d 1 ' 


c dx 

£ 2 a* 
K j 

ay 1 *i> « 



a_ 

dy 



id n 

d_ Pd_ 

fjL.il 


, *0 v 

l** 2 J 

dy p etc 

U**J 

* 2 J 






* 0 . 


(14) 


To deter min e the modal representation of the electromagnetic field we need to solve (14) 
subject to the pseudo-periodic boundary conditions [13] 

rsn H t (d/2,y) - e** d HJ-d/l.y) , 


E t (d/2,y ) = € m9a E t {-dl2,y), 

— (rf/2,y) - e t,od ^(-dl2,y), ^(<f/2,y) » ^(-<f/2,y) . 

etc ax ax ax 


( 15 ) 


In addition, we also need the interface conditions for the fields and their derivatives at the 
medium discontinuities. However, the burden of mentioning the interface conditions can be 
relieved if we understand Maxwell’s equations in the sense of distribution [13]. 


Equation (14) is a coupled system of equations for E z and H v whose direct solution 
seems to be difficult. However, by the following field decomposition, its solution is simplified. 
In a z- and y-invariant medium, an electromagnetic field is said to be Ex (Hx) if the x- 
component of its electric (magnetic) field vanishes. Let the superscript ( e ) denote the Ex 
field, and ( h ) the Hx field. Then from (9) 


i-£« ♦ 

dx ’ 

I* - o. 

ay 

(16a) 

- 

tk,4-E? -0. 

(16b) 

dx 1 




In Appendix A, we prove that any field (E z ,H z f, where superscript T signifies matrix 
transpose, satisfying (14) and (15) can be uniquely expressed as a sum of an Ex field and 
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an Hj. field, and these two component fields satisfy (14) and (15) independently. Thus the 
task of solving the vector-valued boundary-value problem (14) plus (15) is reduced to two 
independent scalar ones to be derived below. 


Using (16a) to eliminate H z from the first equation of (14), and using (16b) to 
eliminate E z in the second of (14), we obtain. 


“1 

( 1 d _(«>' 

| * -£-£?> * Pe? ■ o. 

1 dy 1 

(17a) 

•u 

€ ebr * ) 

* iis* * Ph? = o. 
ay’ • 

(17b) 


In these two equations, the electric field and the magnetic field are no longer coupled. 
Actually, these two equations are identical to the equations for the TE and TM polarized 
fields in classical mountings [5-8], provided that P is replaced by P . 


The symmetry exhibited by (17a) and (17b) with respect to and and with 
respect to e and n suggests that they can be rewritten as 

+ Pf® = o. ( 18 ) 



1 9 pa 

+ 


dx 

o«ax J 


dy 2> 

* e, h, and 



w 

9 9 

F W , ffW 



Let denote quantities complementary to F® such that 

G w = H z \ G (k) = E z h) 

Then, (16a) and (16b) become 


(19) 


(20) 
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(21) 


_i G w = — i — F w 

dy o w *o & 

where 

5 W * - 1 , d w ■ 1 . ( 22 ^ 

The new notations introduced in (19), (20), and (22) allow the Ex and Hx fields to be 
treated identically in the rest of this paper. 




Let a trial solution of (18) be 

F w (x,y) » 

Then the standard procedure of separation 

t 0)) 


i d 

1 d 

dx 

k o w (tc 


««(x) u«(y) . 
of variables leads to 

+ <* 2 -* W2 )« w =0, 


(23) 


(24) 


where k^ is a constant. The differential equation (24) and the boundary conditions (15) 
together pose a boundary-value problem which is considered in the next section. 


V 


2.3. The Boundary-Value Problem 

In this section and the next section, for simplicity, we omit the superscript (s) in the 
relevant quantities. Let L be a differential operator defined by 

as) 

dx\ a dx ) 

Then the boundary-value problem for determining the eigenvalues and the eigen-functions 
of the modal fields is given by 
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(26) 


Lu - p u, 

u (d/2) » e ,m *'u(-d/ 2), • 

u'(dl 2) = e'**V(-<f/2), 

where a prime indicates the differentiation with respect to x, and 

p-X 2 . 

Let us now define an inner product ( , ) for any two integrable, complex-valued functions 
u(x) and v(x) by 

(U.v) = / u(x) v(x) dx , (28) 

J-rf/2 o(X> 

where a bar indicates the complex conjugate. For a lossless dielectric grating (e(x) > 0, ji(x) 

> 0), it is easily seen that L is self-adjoint, i.e., 

(Lu,v) ■ (u,Lv). (29) 

From the theory of ordinary differential equations [14], we know that the eigenvalues 
determined by the boundary-value problem (26) are real and they form a denumerable 
sequence. Furthermore, the eigen-functions form a complete, orthonormal basis in the sense 
that any continuous and piecewise differentiable function /(x) satisfying the pseudo-periodic 
boundary conditions in (26) can be expanded in the eigen-functions. 


In order to embrace the most general cases, however, we assume that region 0 is 

• — * * ■ .fr y - lfry .-- -* -1 - ... 

composed of two media of complex permittivities and permeabilities, i.e., the functions e(x) 
and n(x) are in general complex valued. In addition, we assume and possibly k p to be 
complex. (This is a minor generalization of the works of Ref. 5-9.) This permits us to apply 
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the current model to the homogeneous problem of grating diffraction [15]. In either case, 
the operator L is no longer self-adjoint; therefore, the eigenvalues of (26) are no longer 
necessarily real and the eigen-functions are no longer orthogonal and complete. In order 
to be able to use the modal field expansions for the total field, it is necessary to consider the 
adjoint of (26), which is defined by 

L*u * * p*u\ 

H*(rf/2) * ■ (30) 

u-'(d/2) - t r '- i u-'(-dl2), 

where the superscript + indicates the adjoint and L + is the differential operator adjoint to 
L. It is easily seen that 

(I«,V) « (u,L*v*) < 31) 


if 



( 32 ) 


From the theory of non-self-adjoint boundary-value problems [14,16], we know that under 
certain conditions, which (26) and (30) satisfy, two mutually adjoint boundary-value problems 
have the following properties: (a) Both boundary-value problems possess an infinite number 
of eigenvalues and the eigenvalues can be ordered such that 

P ; - P„. m-0,1,2..... < 33) 


(b) The eigen-functions {u m } and {««}, are bi-orthonormal, i.e., 


( 34 ) 


(c) Any continuous and piecewise differentiable function f(x ) satisfying the boundary 


conditions has a uniformly convergent formal expansion 

fix) - £ </,0 «„(*)■ ( 35 ) 

m -0 

Hence, even for lossy dielectric or metallic gratings, it is still mathematically justified to 
represent the total electromagnetic field in region 0 by a superposition of modal fields, as 
has been done by Botten et al. 

Incidently, the Rayleigh expansions (10a,b) and (lla,b) can be viewed as expansions 
in basis functions 

*„(x) - ( 36 ) 

It is easy to verify that e n (x) are eigen-functions of (26) with L replaced by cP/dx 2 . The 
adjoint of this new boundary-value problem, with respect to a new inner product ( , ) 
defined by 

(«,v> - f*\(x)v(x)dx, . . ( 37 ) 

is (30) with L + again replaced by (P/dx 2 , and 

«;w - ^" v ( 3g ) 

a 

are the adjoint eigen-functions. 

2.4. Eigenvalues and Eigen-Functions 

The explicit forms of the characteristic equation for determining the eigenvalues and 
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eigen-functions can be most conveniently derived by taking advantage of the simplicity and 
symmetry of «(*) and n(x) given in (1). It is easy to verify that the following two functions 
are two linearly independent solutions of (24) 




COSY!*, 

Y,4 


COS- 


COSY 2 


1 * 1 -!) 


°* Y » ci„ Y ^««v Ll- 
sin-— sin y 2 |*| - — 

°iY 2 2 V 2 ) 


Os|*|s^, 


7‘ w ‘f 


(39) 


— dnY!*, 

Yi 

— sgn(x) 
Yi 


• Y l^i 

sin cos y ? 

2 


. dA o 2 Yj Yi^i . 

*| — ^ cos— J-^smYj 
2 J ffiY2 2 


w-4 

1 1 2 


Os|*l*y. 
d x , - d 
(40) 


where, for j = 1, 2, 

2 rl (41) 

Yy * fy - P • v 

Clearly, <p e is an even function and <p 0 is an odd function. The general solution of (24) is 

therefore given by 

«(*) -i4<p.(*) + B*,(x), ( 42 ) 

where A and B are arbitrary constants. Imposition of the boundary conditions in (26) on 
(42) gives the characteristic equation 


A(p) 


( 1 - k) <p # (d/2) ( 1 + k) y 0 (dl 2) 

( 1 + k) q>'(<//2) ( 1 - k) «>'(<// 2) 


» 0, 


(43) 


where 


It is obvious that when k = ± 1, i.e., in Littrow mountings for which the normal incidence 
is a special case, A(p) is a product of two factors. For general angles of incidence, 
factorization of A (p) is impossible, and the explicit form of the characteristic equation is 


cosy^cosy^ - + — ^ ) sinY sin y^ - cosa^ = 0. 


2 i°iY 2 Wi ) 


(45) 


This is a transcendental equation for p, whose solutions are, in general, complex numbers. 


The eigen-functions of (26), expressed in terms of <? e and <p 0 , are given by 

iC<p fl ( x), tf k* + 1. <P„(<*/2) =0, 

«(x) * »C<p 4 (x), if 1, *,(<*/ 2)=0, 

C [ ( 1 +k) <p,(d/2) <p,(x) - ( 1 - k) 9 ,<d/ 2 ) *,(*) ] , else, 

(46) 


where C is the normalization constant. The eigen-functions of the adjoint problem (30) can 
be simply obtained by replacing the relevant quantities in (39-46) by their adjoint 
counterparts. It can be shown that for each eigenvalue, there is in general only one eigen- 
function and accidental degeneracy of an eigenvalue can only occur in Littrow mountings. 
Since the normalization constants C for u(x) and C + for u + (x) are not individually fixed by 
(34), we can demand = C. Then it can be shown that 

u*(x) = k _ 1 «(-x). 


This direct relationship between the two mutually adjoint eigen-functions is very useful in 
the numerical implementation of the theory. 
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2.5. Modal Field Representation 

Since the basis functions {u m w } are complete and bi-orthogonal in the sense stated 
in 2.3, the general solution of (18) for F® can be written as 


where 


F 0) = £ »?(>)«?(*), 

*-0 


<•>«(?) 3 + 


(48) 


(49) 


and aj 3) and b m (t) are modal field amplitudes to be determined later. Substituting (48) into 
(21) and integrating with respect to y, we obtain an expression for Gfi\ 


- £ Xm(y)"m( x )' 

M‘0 


where 


x?00 * 


k: * x 


W2 


,(*) 


-[ - b®cos\ ( £y + a? sinXjfy ] 


(50) 


(51) 


and 


w“(l) = 


1 *« s® 




(52) 


(The integration constant that would appear in (50) can be shown to be zero.) Thus, by 
(48), (50), and (A-l) we complete the derivation of the following modal field representation 
of the total electromagnetic field in region 0: 


E t (x,y) ' 
[H t (x,y) J 


■ E 

M *0 




x?(y )*?(*) 


(53) 
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Before closing this section, we give the orthogonality relation between the vector- 
valued Ex solution and Hx solution. For this purpose, we define an inner product [ , ] of 
two vector-valued functions 9 and ^ such that 


rW 

a / <P r 

J - dll 


f-i/i* 0 

0 ut) 


1 ) 1 dx. 


(54) 


Let = (u m ^, w m W) T , and u+^) T . Then it can be shown that 


f -(«) _ t W'| A 

IV**9 b J ■ °- 


(55) 


2.6 Matching Interface Conditions 

Having obtained the expressions of the electromagnetic fields in all regions of the 
space, we are now ready to form the final system of linear equations for determining the 
unknown field amplitudes by applying the interface conditions aty = ± h/2. The interface 
conditions are the following: 


E = E , 

1* 1- ’ 

H z . = H t _, 

(56a) 

E » E , 

x* x- > 

- H x _, 

(56b) 


where the subscripts ± indicate limits from above and below the interface respectively, and 


from (9a) and (9b), 


E x - — I 

'k±E t « 
k Z cbc 1 

■ 

(57a) 

H x * — I 

k—H 

- .k,±E t ). 

(57b) 

* £ 2 ' 

l a* 1 



Substituting (10a,b), (lla,b), and (53) into (56a) and (56b) and carrying out some tedious 
algebra, we have, for the continuity of E z and H z aty = + h/2. 
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(58a) 




‘*w ♦ E *“«.(*> 


and 


e [(«.*’♦ Ww * w ^rcK’w 


m -0 

?(*) 


(58b) 


r«o(^) ♦ e *?’ ».(*)• 


For the continuity of £ z and // z at y * - /t/2, we have 


and 


E*?‘.w < 58c) 

m-0 


i[(tf-wo‘W*ww]- < 58d > 

m "0 m - -• 


For the continuity of £ x and H x at y = + Zi/2, we have 

jw-0 ^ 


(59a) 


= ( tf> p? i?> ) S W - E ( «. ) «.(*> 


and 


m-0 € 

K’C/f - t?sD-«iC) - £ 


(59b) 


Finally, for the continuity of E x and H x at y - - h/2, we have 
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and 


- E * »?«.*?)«.<*) ( 59c > 

* -0 ^ »■-• 




««>(*) 


E (*?»?*? 


<?«.?)**>• 


(59d) 


In the above equations the unknowns are 


- 0 «co.(*2>A/2), 


(60) 

i . Rf e ‘ t '" tn , 

If - Tfe"> 1 , 

(61) 


and the rest of the new notations are defined in Appendix B. The numerical solution of 
equations (58) and (59) will be considered in section 3.2. 


2.7 Diffraction Efficiency and Polarization 

Once Rj*) and are solved from (58) and (59), Rj- S) and T n {s) are given by (61). 
By virtue of (9a) and (9b), all quantities of practical interest can be readily expressed in 
terms of R n ^ and Suppose the media in regions 1 and 2 are lossless, and the incident 
plane wave is normalized such that 



(62) 


Then, the diffraction efficiencies for the reflected and the transmitted propagating waves of 
order n are given by 
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(63a) 


ft (i) 

a) P* 




and 


If the media in region 0 are also lossless, the energy balance theorem holds: 


0) 


(63b) 




(64) 


where n and n ' run through all propagating orders in regions 1 and 2, respectively. 


In many applications involving conical mountings, it is very important to be able to 
predict the states of polarization of the diffracted orders. Let us associate with a 
propagating order having a wave-vector 

± + zk t , 


two unit vectors $ n ® and p® such that 


i, 


( 7 ) 


sh i 


r</) 


ra * frU) 


< 7 > .. 


( 66 ) 


In (65) the plus sign is for j = 1, and the minus sign is for j = 2. The polarization of this 
diffracted order can be described with the following two angular parameters [17]: 
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v 



2 


. 0 ) 


tU) 


arctan 


\K 


<j) i 


l£ a) 

i ) 


-arg 


' vW 

L nt 

pO) 

V £ v / 


- it < <£ it , 


( 67 ) 


O') 


where £„/) an d £^0) are the 5- and p-components of the electric field. (The notation a n 
should not be confused with a„ defined in (12).) It is a simple exercise to show that for a 
diffracted wave in medium 1, 

-(l)V dW . od)t D® 


.0) 




(68) 


C N € (1) 

The expression for a diffracted wave in medium 2 may be obtained from (68) by replacing 
£„ (,) by r„ (,) , superscript (1) by (2), and by - 0„ (2) . 


3. NUMERICAL ASPECTS 

3.1 Solution of the Characteristic Equation 

An efficient, reliable, and accurate numerical method for solving the characteristic 
equation (45) is of vital importance to the performance of a computer program 
implementing the modal analysis. Two very different numerical methods have been 
employed. The method of Botten et al [7,18] is a general one capable of finding all zeros 
of an analytic function in a prescribed region of the complex plane. That of Suratteau et al. 
[8] and Tayeb and Petit [9] is a problem-specific method that takes advantage of the fact 
that (45) can be factored in Littrow mountings. Both of these methods systematically find 
all eigenvalues of (45) in a prescribed region of the complex plane, and both of them 
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perform well even for highly conducting gratings. As reported by the authors of the second 
method, their method is as effective as the first method but requires significantly shorter 
computation time. For this reason, we have decided to use the second method in our 
numerical implementation of this work. For details of these numerical methods, the 
interested reader is referred to the references cited above. 

3.2 Solution of the Field Amplitudes 

Before embarking on the numerical solution of (58) and (59), we first analyze the 
composition of these equations. In (58) and (59) there are eight equations and eight sets of 
unknowns. The right-hand sides of these equations are expanded in basis functions {e n }, 
which are bi-orthogonal to their adjoint {e*}. The left-hand sides are expanded in four 
different sets of functions {u m ^}, {u m ^}, {w m ^}, and {w m ^}. Of these four, as scalar- 
valued functions, {u m (e) } and {uj h) } are bi-orthogonal to their respective adjoints, but not 
to each other. The functions {w m (<) } and {wj h) } are proportional to the derivatives of 
{uj e) } and {uj h) } and they are not orthogonal to any other functions. 

Equations (58) and (59) constitute a system of equations in known function 
expansions with unknown expansion coefficients. Such a system can be solved by the 
method of moments [19], which consists of three steps. First, a projection basis, i.e., a set 
of linearly independent testing functions, is chosen. Then, both sides of the series expansion 
equations are projected onto this basis by forming appropriate inner products with the 
testing functions. This step eliminates the x-dependence of the equations and produces an 



algebraic linear system of equations of infinite dimension. Finally, the linear system is 
truncated to a finite order and its solution is obtained by the standard numerical techniques. 

Since equations (58) and (59) are already expanded in terms of bi-orthonormal basis 
functions, it is advantageous to choose the adjoints of these basis functions as the testing 
functions so that the subsequent numerical solution can be simplified. For each of the eight 
equations we have two convenient projection bases, and {e + }. Therefore, there can 

be many different combinations of choices of projection bases for the overall system. 
Following Suratteau et al , a projection method in which the interface conditions for the z- 
components are projected onto one basis and those for the x-components are projected onto 
the other basis is called a hybrid method. A projection method in which all interface 
conditions are projected onto one basis is called a homogeneous method. For the case of 
non-conical mountings, Suratteau et al have proved that the numerical solutions resulting 
from the hybrid methods satisfy the energy balance and reciprocity criteria automatically 
(independent of the truncation orders) while those resulting from the homogeneous methods 
do not. It can be shown that the above statement is also true in the case of conical 
mountings. 

In this paper, we adopt the homogeneous method utilizing the projection basis 
« w >. so that the energy balance and reciprocity criteria are not automatically satisfied. 
Multiplying equations (58) by d^ s \x)/a^\x), and equations (59) by Q^ s \x), then integrating 
over a grating period, and making use of the bi-orthogonality relation (34), we have the 
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V 


V 


V 


V 


w 


following matrix equations: 


WY - UX + l//, 
Dr * QX + P/, 


(69) 

(70) 


where 


r ' 


' ' 


fi*) * ) 

l z 5 0 - 

fL*) 

, Y * 


, / * 

0 

a* 


a<*> 

" ) 


'T^o. 

, 0 ; 


5^ is the Kronecker delta, and the rest of the matrices are defined in Appendix B. In (7 1 ), 


each element of the column vectors is itself a column vector, and the elements of A!" and Y 


are related to the unknown field amplitudes listed in (60) and (61). The matrix D in (70) 
is diagonal, so vector Y can be expressed in terms of X without numerical matrix inversion. 

Substituting the expression of Y into (69), we have 

(WD~'Q - U)X * (U - WD' x P)l. ( 12 ' ) 

This is the final linear system of equations from which we numerically determine the field 


amplitudes. 


In order to solve the linear system (72) on a computer, we unavoidably have to 
truncate the matrices. We designate N as the total number of terms retained in Rayleigh 
expansions (we truncate the Rayleigh expansions symmetrically with respect to the zero 
diffraction order) and M as that retained in the modal expansions. The integers N and M 
are called the truncation orders. It is easily seen that for the solution of (72) to be well 
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(neither under- nor over-) specified, the two truncation orders must be the same. Thus, the 
matrices W y U, P, Q, and D in (72) are 4N by 4 N square matrices. 

4. NUMERICAL EXAMPLES 

In this section we present some numerical results. The computer program is written 
in Fortran 77 and double precision is used for real and complex arithmetic. For the special 
case of non-conical mountings, the program has been checked using published data with 
good agreement. For the general case, it meets the energy balance and reciprocity criteria 
with reasonable accuracy. Table 1 tabulates TE and TM diffraction efficiencies of a metallic 
gra ting in a non-conical mounting. The data of Botten et al are taken from Table 1 of Ref. 
6. Note that Botten et al adopted a hybrid projection method that allows unequal 
truncation orders N and M. Clearly, the agreement is very good, especially for the TE 
polarization. 

To date, there is no numerical data, especially data of the polarization parameters, 
for non-perfectly conducting lamellar gratings in conical mountings available in the literature; 
therefore, we present some original data in the rest of this section. Listed in Table 2 are 
diffraction efficiencies (rj), polarization angles (a, 6), and diffraction angles (0,0) of a 
dielectric grating in a conical mounting (the parameters are listed in the table caption). 
Note that the incident plane wave is right-hand, circularly polarized. Listed in Table 3 are 
diffraction efficiencies (rj), polarization angles (a, 6), and diffraction angles (0,0) of a highly 
conducting grating in a conical mounting. The incident plane wave is linearly polarized with 
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equal s and p component amplitudes. In both cases, the truncation orders are chosen to 
ensure that the accuracy of the data is better than one percent (see discussion below). 


In Figs. 3a, 3b and 3c we show the change of diffraction efficiencies rj and 
polarization angles a and S of a small period dielectric grating as the incident azimuthal 
angle <p sweeps through the first quadrant. The incident plane wave is always p polarized 
and it strikes the grating from the optically denser medium at a polar angle greater than the 
critical angle (for total internal reflection). This configuration is reminiscent of what occurs 
in a planar waveguide grating coupler for a TM polarized guided-wave [20]. The reflected 
and the transmitted negative first orders pass off at about <p = 37° and <p — 60°, 
respectively. It is evident that as soon as 0 is nonzero, the diffraction orders become 
elliptically polarized. The two first orders are nearly circularly polarized at 0 - 20°. Also, 
the senses of polarization of the two first orders remain right-handed throughout the angular 
range of their existence. 


Next, we illustrate the excellent convergence rate of the modal method. For this 
purpose, we define a measure of error as follows 


A* = log 


to 


f*-r 

r 


(73) 


where f N stands for any one of the physical quantities, such as a diffraction efficiency or a 
polarization angle, computed with truncation order N, and / = / N ., where N* > N is an 
integer. If/ stands for the sum of the diffraction efficiencies for a lossless grating,/ = 1.0. 



Roughly speaking, the negative of gives the number of correct significant digits in the 
numerical results. 

Fig. 4 shows the convergence of the sum of the diffraction efficiencies and the 
efficiencies of the negative first orders in reflection and in transmission for a dielectric 
grating. The physical parameters are the same as those of Table 2. The truncation order 
N varies from 11 to 61 in increment of 2, and N* - 63. It is evident from the figure that 
better than one percent accuracy is achieved as soon as the truncation order is greater than 
11. The convergence is not monotonic. The large oscillation in the convergence sequence 
of the reflected order is probably due to the smallness of the diffraction efficiency (see row 
3, column 2 of Table 2). If we make a low order polynomial fit of each set of the data in 
the figure, the three resulting curves will have more or less similar shapes and close 
locations. This implies that the energy balance criterion can be used as a good accuracy 
indicator, thanks to our choice of the homogeneous projection method. 

Fig. 5 shows the convergence of the diffraction efficiency and diffraction angles a and 
6 of the negative first order of a metallic grating. The physical parameters are the same as 
those of Table 3. The truncation order N varies from 11 to 69 in increment of 2, and AT = 
71. Since now the grating is metallic, the convergence in this case is, as expected, slower 
than that of Fig. 4. However, better than one percent accuracy can still be achieved with 
a truncation orders of 40 or greater. 
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As the authors of Ref. 5-8 have shown, one of the unique features of the modal 
method is that it works very well even for deep, highly conducting gratings. This feature is 
also true for the extension of the method to conical mountings as developed in this paper. 
This is understandable, since the fundamental analytical and numerical issues for conical and 
non-conical mountings are essentially the same. For the sake of saving space, however, we 
will not provide any numerical evidence here. 

5. SUMMARY 

In this paper, we have extended the rigorous modal method of Botten et al to the 
case of conical mountings. A crucial step in accomplishing the extension is the field 
decomposition discussed in section 2.2 and Appendix A. The field decomposition reduces 
the vector-valued boundary-value problem given by (14) and (15) to a scalar one given by 
(26), thus tremendously simplifying the subsequent analysis and allowing the previous works 
of Botten et al and Suratteau et aL be used here. The completeness and orthogonality of 
the modal fields in the corrugated region are carefully established. The computer program 
implementing the extended modal method can treat a plane wave of arbitrary angle of 
incidence and polarization. It converges very well for highly conducting grating materials 
and very deep gra tin g grooves. We have included some original numerical data of both 
diffraction efficiencies and diffraction polarizations for conical diffraction configurations. 

The mathematical formulation presented in this paper has been kept general. In fact, 
the explicit x-dependence of c(x) and n(x) given in (1) is not used, except for the derivation 



of the characteristic equation and the eigen-functions in section 2.4. Therefore, the 
formulation is valid for any periodic e(x) and p(x). However, it is precisely the simple x- 
dependence of e(x) and p(x) in (1) that makes the easy derivations in section 2.4 possible. 
For any other permittivity and permeability variations, the solutions for the eigenvalues and 
eigen-functions become very complicated, and the modal analysis quickly loses its advantages 
over other grating methods. 

The present analysis can be easily extended to treat gratings of arbitrary groove 
shapes and waveguide gratings with a number of uniform layers above and below the 
corrugated region. These extensions will be the subject of a future paper. 
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7. APPENDICES 

7.1 Appendix A 

Theorem: 

Any solution of (14) and (15) can be decomposed such that 




' E®' 


E™' 


s 


+ 


"J 




u (k) 

1"* J 


where (E z (e) , H z (e) ) T and (E z {h) , H z (h) ) T are Ex and Hx, respectively, and they satisfy (14) 


262 



and (15) independently. The decomposition (A-l) is, in general, unique. 
Proof: 


Since there is no explicit ^-dependence in (14), any solution of (14) is necessarily of 
the following form 




' E^ix.y) ' 
k E zx( x »y) ; 


(A-2) 


where 


r E^x,y) ' 
k H^y) 


= g “y 


' **(*) ' 


(A-3) 


Substituting (A-3) into (14), and eliminating the ^-dependence, we have 


£l±(-L± 

« dx\?dx) 




a 1 - x 1 


K P 

+ iX-i — 
*b « 


(— — 1 






+ * 2 - X 2 

dr * 2 A) 






= o. ( A - 4 ) 


Equation (A-4) is a linear, two-dimensional vector-valued, second-order ordinary differential 
equation. It, therefore, has four linearly independent solutions. Suppose [E zV H zX ) J is E. 
Then, from (16a) 

(A5> 


Substitution of (A-5) into (A-4) results in two second-order equations for l x (x). These two 
equations are not independent; one can be derived from the other. Thus, we have two 
linearly independent solutions of (A-4) that satisfy the Ex condition (16a). Similarly, if we 
demand (£ zi , H zX ) T to be Hx, then 



iXekqdx 


(A-6) 


f A C ^ (M 

and we have two linearly independent solutions of (A-4) that satisfy the Hx condition (16b). 
Gearly, the solutions of Ex type and Hx type are linearly independent. Therefore, any 
solution of (A-4) has the following form 


' i,w ' 

(*) 

1 

"XI 

^ (<) 

U X2 

, -<*) 

vg ] 

^ „(A) 

v » ) 
v 12 

> *•»<*> , 

* Cxi 

vff : 

v *1 ; 

+ 

v « 
V « J 

+ C At 

«g j 

+ C A2 

u (h) 

“u ; 


(A-7) 


V ) \ V£, ) y va ) yu n ) y u u ) 

where c v ^\ l = 1, 2, s = e, h, are constants, (uj- e \ v u w ) T and Uj} h) ) T are Ex and Hx 
solutions of (A-4) respectively. So, the decomposition (A-l) is always possible. 

Next, we impose the pseudo-periodic boundary conditions (15) on the general 
solution (A-7). This leads to the following characteristic equation 


u xi 

TjV 

U X2 

Kg 

Kg 

u { ‘ y 

u x\ 


u®/ 

Kg' 

iM 

Y x\ 

Kg 

r/W 

'■'Xl 

ug 

Y Xl 

lArf 

Y X1 

^Xl 

tfg 


= 0, 


(A-8) 


where 


l/g - «®(<//2) 

Kg ■ »{J(<//2)-«'*^v‘J(-J/2), 

t/g' - «?/(<*/2) #g(-<f/2). 

Kg' - v«'(d/2) v“»'(-<</2). 


(A-9) 


By elementary row manipulation and making use of (A-5) and (A-6) for the Ex and Hx 
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solutions, it can be shown that the two off-diagonal two-by-two matrices in (A-8) can be 
made zero. Therefore, the Ex and Hx solutions of (A-4), and hence those of (14), satisfy 
the boundary conditions (15) independently. 

Suppose the decomposition (A-l) is not unique. Then we may have a decomposition 
of the zero field into two non-zero orthogonal modal fields. Furthermore, each of these 
fields satisfies both the Ex and the Hx conditions simultaneously. However, this leads to 
k z 2 + X 2 = 0, a condition which is, in general, not true. This completes the proof of the 
theorem. 


7.2 Appendix B 

As in the main text of the paper, j = 1, 2, and s = e, h. In (58) and (59) 

, j(f) g -itfw (B-l) 


w * 2 ♦ *? 2 X® h 
A? = — — tan-^- 

■ i<!) 2 


m . * 2 + x; )2 _,_ o 


B? = 




ctan- 


(B-2) 


.CO 




2 € 


0) 


,U) , 




(i) 


0)2 




a l — — . 

* k u)2 


(B-3) 


tun • ? 

For the sake of clarity, the matrices in (69) and (70) are expressed in block forms, 
followed by the definitions of each sub-matrix. 




W = 


1 

1 

W< 4) A (k) 

W il) B (h) 

1 

- 1 

- W^A® 

jy«) £<*) 


w™ b m 

1 

1 

_ ^a™ 

jy(h) £(*) 

1 

- 1 


(B-4) 
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k 0 0 
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o (t) 

0 

0 

0 

' -T?0® 


r (t) « 
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42-4?«.. 


(B-8) 



(B-9) 



(B-10) 


- (W« O,. 

(B-ll) 


To get the expression for matrix P, we only need to reverse the sign of in (B-7). In (B- 
10) and (B-ll) the inner products are those defined in (28) and (37), and the subscript s 
indicates that the weight function a (1) should be used in the integral. By using (55) and (47), 
it can be shown that the two sub-matrices in (B-ll) are related such that 

+ *£ = °- ( B ‘ 12 ) 
Hence, there are five sub-matrices involving the eigen-functions to be calculated. It can be 
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shown that these sub-matrix elements can be expressed in terms of the left and right limits 
of u m (J) and its derivative at ± d x l 2, if the periodic medium is characterized by (1). 

Since and r 3 ^ are proportional to k p the matrices defined in (B-4) through (B- 
7) become block-diagonalized when k z = 0. Of course, this means that in non-conical 
mountings, E z and Hj are de-coupled. 
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TABLE CAPTIONS 


1. Numerical comparison with the results of Botten et al for a metallic grating in a non- 
conical mounting. Parameters: d - 1.0 /xm, d x = 0.4001 nm, h - 0.1 ^m, e (1) = c (2) 
= €j = 1.0, c 2 = (1.5+/1.0) 2 , A, 0 = 0.8 jxm, 0 = 11.5°, 0 = 0°. Truncation orders: 
Botten et al, M = 20, N = 51; this paper, M = N — 31. 

2. Diffraction efficiencies (rj), polarization angles (a, 5), and diffraction angles (0,0) of 
a dielectric grating in a conical mounting. All angular values are in degrees. 
Parameters: d - 1.0 /xm, d x = 0.5 /urn, h = 0.5 jxm, e (1) = e x = 1.0, e (2) = e 2 = 2 . 25 , 
A. 0 = 0.5 nm. Incident polarization: a = 45°, 5 = 90°. Incident angle: 0 = 0 = 45 °. 
Truncation orders: M = N = 31. 

3. Diffraction efficiencies (rj), polarization angles (a, 6), and diffraction angles (0,0) of 
a metallic grating in a conical mounting. All angular values are in degrees. 
Parameters: d = 1.0 /xm, d x = 0.5 /xm, h = 1.0 /xm, € (1) = e x = 1.0, e (2) = e 2 = 
(0.1+/5.0) 2 , A. 0 = 0.5 /xm. Incident polarization: a = 45°, S = 0°. Incident angle: 
0 = 30°, 0 = 45°. Truncation orders: M = N = 51. 
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FIGURE CAPTIONS 


1. The coordinate system for a lamellar grating in a conical mounting. 

2. The geometry of a lamellar grating. 

3. Diffraction efficiencies rj (Fig. a) and polarization angles a (Fig. b) and 6 (Fig. c) of 
a dielectric grating versus the incident azimuthal angle 0. The polarization angles are 
measured in degrees. Parameters: d = 0.3 /mi, d l = 0.15 /mi, h = 0.15 /im, e (I) = 

= 1.5, c (2) = e 2 = 1.0, X 0 = 0.5 /tm. Incident polarization: a = 0°. Incident polar 
angle: 6 = 60°. Truncation orders: M - N = 45. 

4. Illustration of convergence for a dielectric grating. Round dots: sum of all diffraction 
efficiencies. Hollow squares: diffraction efficiency of - 1 order in reflection. Hollow 
triangles: diffraction efficiency of - 1 order in transmission. The parameters are the 
same as for Table 2. 

5. Illustration of convergence for a metallic grating. Round dots: diffraction efficiency 
of -1 order. Hollow squares: polarization angle a. Hollow triangles: polarization 
angle S. The parameters are the same - as for Table 3. 
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Fig. 3c 
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